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ANALYSIS OF CARBON NANOTUBE PULL-OUT FROM A POLYMER MATRIX* 

S. J. V. Franklandand V. M. Harik 


Abstract. Molecular dynamics (MD) simulations of carbon nanotube (NT) pull-out from a polymer 
matrix are carried out. As the NT pull-out develops in the simulation, variations in the displacement and 
velocities of the NT are monitored. The existence of a carbon-ring-based period in NT sliding during 
pull-out is identified. Linear trends in the NT velocity-force relation are observed and used to estimate an 
effective viscosity coefficient for interfacial sliding at the NT/polymer interface. As a result, the entire 
process of NT pull-out is characterized by an interfacial friction model that is based on a critical pull-out 
force, and an analog of Newton’s friction law used to describe the NT/polymer interfacial sliding. 

Key words, carbon nanotube, polymer nanocomposite, interface 

Subject classification. Computational Materials Science 

1. Introduction. Polymer nanocomposites with carbon nanotubes (NTs) as fillers have been 
fabricated and mechanically tested to explore their potential as strong, lightweight materials [1, 2], While 
some reinforcement of the polymer matrix has been observed in the form of enhanced moduli relative to 
the polymer matrix alone, a detailed understanding of the NT/polymer interfacial region remains 
fundamental to the optimization of the mechanical behavior of these composites. The interface plays a 
significant role in the stress transfer and the consequent improvements in the nanocomposite stiffness and 
strength. Characterization of the NT/polvmer interfacial interactions during pull-out will be the focus of 
this study. 

For traditional fiber-reinforced polymer matrix composites, the interfacial shear strength is typically 
evaluated by employing fiber pull-out [3] or fiber push -out tests [4], The fiber decohesion process is 
characterized by the critical shear stress required to debond the fiber. The value of the shear stress is 
estimated from the critical axial load using the shear-lag model. For nanocomposites, a similar procedure 
can be utilized both experimentally [5] and theoretically [6], For a NT embedded in a polymer matrix, the 
critical shear stress at the interface has been evaluated via molecular dynamics (MD) simulations and an 
analog of the shear-lag model [6], A comparison was made between nanocomposites with and without 
the nanotube chemically bonded to the polymer indicating that chemical bonds reinforce the interface [6], 
This study extends the previous work for the case of a non-bonded NT to refine the analysis of interfacial 
sliding between the NT and polymer nanostructured surfaces. 

In the system simulated here, the matrix is crystalline polyethylene, which interacts with the NT via 
van der Waals forces represented by the Lennard-Jones potential. MD simulations are used to confirm the 
critical pull-out force and investigate interfacial sliding as a NT is pulled through the surrounding 
polymer. NT velocity and displacement are traced in order to characterize the interfacial interactions 
during NT pull-out. The NT/polymer interfacial friction is described with a relation between the applied 
force to the NT and the NT velocity that is an analog of the Newton’s viscous friction law, which involves 
an effective viscosity coefficient of the interfacial sliding. Furthermore, an interfacial friction model for 
the entire pull-out process is proposed. 
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2. Molecular Structure of a Nanotube/Polymer System. The nano-scale material system in this 
study is composed of a carbon NT and a polyethylene (PE) matrix. Specifically, a (10,10) carbon NT is 
embedded in crystalline polymer matrix (Figures 1 (a) and (b)). The system contains about 23,000 atoms 
including an 880-C-atom NT and 178 chains of 42 methylene units with the hydrogen atoms explicitly 
included. The dimensions of the system are ~6 x 6 x 6 nm. In the molecular dynamics (MD) simulations, 
the entire system is represented with the hydrocarbon potential developed by Brenner, et al. [7], The 
system is replicated across periodic boundaries in each of the 3 -dimensions, thereby making the NT and 
the PE chains infinitely long. This composite is similar to the crystalline non-bonded polyethylene one 
used in the previous work [6], and the detailed chemical description is kept for consistency. 



FIGURE 1 . Molecular structure of crystal line polyethylene NT composite: (a) a cross-section of centra i 3. 6 nm of the 
sample (b) a side view of NT and two nearby polyethylene chains. 


The entire NT pull-out process is simulated by MD. As the sample cell is periodic, the pull-out is 
simulated by pulling the NT through, rather than out of, the PE matrix. A unidirectional force is applied 
to each atom of the NT along the NT axis. The pull-out force is increased incrementally over time 
(Figure 1). The simulation is run for ~ 100,000 MD time steps of 0.5 fs. Data generated from MD 
simulation of this system forms the input to the interfacial friction model described below. 

3. Characterization of Nanotube Pull-out and Numerical Results. The MD simulation captures 
the full NT pull-out process. Initially, after a sufficient amount of force is applied to the NT in the 
simulation, interfacial sliding begins. In Figure 2, the increments of applied force are shown, and over 
time, the velocity (Figure 3) and the displacement of the NT are monitored (Figure 4). Following the 
observed variations in the axial velocity of the NT (Figure 3), the entire pull-out process can be divided 
into several stages. In the simulation, the axial velocity of the NT is sufficiently greater than the velocity 
of the surrounding polymer in the vicinity of the NT, and the latter is therefore neglected in the analysis. 
A detailed analysis of the NT velocity data is subsequently carried out to characterize the pull-out 
process, in general, and interfacial sliding, in particular. 
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FIGURE 2. The increments of axial force applied to the nanotube over time. 



FIGURE 3. Velocity of the center of mass of a carbon NT during the simulated pull-out. 
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Time steps (0.5fs) 


FIGURE 4. Displacement of individual carbon atoms of the NT during the pull-out process. 

Prior to the onset of sliding, the NT experiences thermal vibrations and interacts with the adjacent 
polymer chains in equilibrium at 300 K. The effect of these interactions on NT axial velocity is 
illustrated in Figure 3 for the first 27,000 time steps (up to time of pull-out, t pu n_ out ) as the applied force is 
being gradually increased to 0.1 nN. Beyond this critical force, the NT starts to slide through the 
polymer matrix. Following the details of MD simulation data presented in Figures 3 and 4 as well as the 
established conceptual framework for viscoplastic interfacial yielding [8], the entire NT pull-out process 
can be described by four stages. 

1) Initially, an equilibrium exists between the NT and the surrounding polymer chains (7 < t pu ii- ou t ~ 
27,000 time steps). In this equilibrium configuration, there exists a net strength of interaction 
between the polymer and the NT which must be overcome for pull-out to occur. 

2) At an applied force of about 0.1 nN, accession of the energy/force barrier occurs initiating 
slippage between the NT and the polymer. This onset of slippage is marked by a noticeable 
increase in the axial velocity of the NT, w hich is observed at t = t pu ii- ou t -27.000 time steps. 

3) Slippage of the NT in the polymer cavity fully develops during a transition stage between t pu ii-out 
and 7 v, v ~5 0.000 time steps, where t ss refers to the onset of steady sliding. 

4) The steady sliding of the NT (for t > t ss ) is characterized by semi-periodic variations in NT 
sliding velocity as the NT accelerates under step-wise increases in the applied force. The onset of 
periodic accelerations in the NT sliding is marked by dots in Figure 3. 

Similar stages of interfacial sliding have been observ ed for other surfaces [9], 


4 



The observed periodic variations in the NT axial velocity (Figure 3) can be related to the size of the 
carbon rings in the NT periodic structure. From a comparison of the velocity and displacement data in 
Figures 3 and 4, a correlation is evident between the distance traveled by the NT and the period of the 
velocity oscillations. The movement of the NT can, therefore, be scaled by the width of a carbon ring, a. 
In Figure 3, the dots on the velocity trace mark the onset of accelerations, which are related to the ascent 
of local energy barriers at the beginning of the carbon ring period. These onsets occur sooner as the NT 
kinetic energy increases relative to the local energy barriers to NT sliding. In the displacement trace 
plotted in Figure 4, there are noticeable changes in the slope marked by the dotted lines, which occur 
simultaneously with the beginning of the accelerations in the velocity trace. These changes in slope are 
separated by 2.46 A, the width of carbon rings, a. 

4. Interface Model and Nano-scale Friction Laws. For fiber-reinforced polymer matrix composites 
undergoing large deformation, interfacial debonding and sliding can be simulated by elastoplastic, viscous 
or viscoplastic models for matrix materials accompanied by a phenomenological interface potential for 
the fiber-matrix interaction [8], A Bingham -type model [8] for viscoplastic failure of the interface relates 
the shear stress, r, to the strain rate, / , via the Newton's friction law (r= uy ) after the yield stress, To, is 
reached: 

(4.1) t = t 0 + yij . 

Here, /u and /u m are viscosities of a fluid and a viscous matrix material, respectively. Newton's friction 
law [10] is defined for the adjacent layers of a continuum fluid that is subjected to shearing. At the 
nanometer level, an analog of this friction law can be applied, which uses the forces applied to the NT and 
the velocity information (Figure 5) to determine an effective viscosity. Furthermore, the stresses and 
strain rates can be recast in terms of average forces and velocity gradients derived from the MD 
simulations of NT sliding. 

In order to characterize and then model the sliding between a NT and the surrounding polymer chains 
on a length-scale level larger than the carbon rings, multiple carbon ring lengths of simulation data are 
needed for sufficiently meaningful averaging of the contributions from local NT structure. The relative 
sliding involves two surfaces, that is, the NT and polymer sliding surfaces. For the longitudinal direction, 
the NT length, F\t, should be sufficiently long, that is, L N j/a >10, where a is the width of a carbon ring. 
The number 10 implies an order of magnitude difference. Likewise, for the perimeter of the NT/polymer 
interface, thc_NT radius should satisfy a similar inequality: 2 7t(R N t+ h vdfv /2)/a >10, where Rm is the NT 
radius and h VdW is an average surface separation caused by the van der Waals forces. The considered 
NT/polymer system satisfies all these averaging (or homogenization [11]) criteria. 

At the length-scale level larger than the carbon rings, an analog of Newton’s friction law (4.1) can be 
developed as follows. The average shear stress, <T rz >, for planar sliding interactions can be defined as 
< T r: > ~ < f > P ui/(Aa), where <f> P uii is the average force applied to the NT and A ss is the common interfacial 
area of steady sliding. The strain rate,/, can be estimated by y~ d<V z >/dr, where <Vy> is the average 
NT velocity in the z-direction and r is the radial distance. Then, the nano-scale analog of Newton's law 
for steady NT/polymer interfacial sliding is 

d < V, > 

(4.2) < >= /V - , 
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where /i e g is ail effective viscosity for the NT/polymer interfacial sliding. As a result, an interfacial 
friction model for the entire NT pull-out process can be developed for the total shear stress, <r rz > pu n, 
which, in addition to the sliding component < T rz >, includes an initial threshold pull-out stress, r 0 : 


(4.3a) 


T rz > pull ~ T 0 + f^eff ' 


d <V„ > 


dr 


After approximating the change in velocity across the NT/polymer interfacial separation, h vdw , as the 
average NT velocity, <V Z > , 


<K> 

(4.3b) < x ri > pull =t 0 + M eff -= » 

"vdW 

where the strain rate is \d<V z >/dr * <V z >/h vdw . The values of the average NT velocity, < V->. and the van 
der Waals separation, h vdw , are detennined from the MD simulation. 

To complete the description of the interface model for the NT pull-out process, the effective viscosity, 
/Jeff, can be detennined from the MD simulation. Such a viscosity coefficient depends on the nature of 
NT/polymer interfacial interactions. In the non-bonded system simulated, the velocity of the NT is, on 
average, linearly related to the average applied force, </>, as can be observed for times longer than t ss in 
Figure 5. The slope of this relation, Xeff, can be parameterized by the MD simulation data. The resulting 
force-velocity dependence 


(4-4) </>= Xeff < V z > 

is an analog of the friction law in equation (4.2). The applied force, <f>, in equation (4.4) is related to the 
sliding component of the shear stress as in equation (4.2) by the force balance: 


(4.5) 


< / >=< A, > A ss 


Therefore, the "effective viscosity coefficient" can be related to the effective viscosity /i e g via 


(4.6) 


Meff ~ 


XeffKdW ! 

2?t{R N T + h v dw ( 2 ) 


Following the analogy, the interfacial friction model for the entire pull-out process is 


(4.7) <f>pull=fo+Zeff <K >> 

where fi> is the critical pull-out force, which is also available directly from the MD simulation. 


5. Characterization of the Interfacial Sliding. The interfacial friction model (4.7) involves the 
critical pull-out fore e,f 0 , and the "effective viscosity coefficient”, ;/«//, the values of which can be obtained 
from the MD simulation of the pull-out process. For every incremental force loading, the NT velocity 
variations shown in Figure 3 are plotted in Figure 5 as a range represented by a vertical line. These ranges 
correspond to NT sliding velocity once the critical pull-out force, To, at 0. 1 nN is reached. After this point, 
the force-velocity dependence can be approximated by a linear relation (the solid line) [9], The inverse of 
the slope of the resulting line is the coefficient Xeff, which is 1.5±0.3 (nN ps)/A. This slope can be 
estimated from dashed lines like those shown in Figure 5. which bound the NT velocity variations during 
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NT sliding. Once the value of the coefficient Xeff, is determined, the effective viscosity, Li e g , can be 
calculated from equation (4.6) as 0.2 cP (centi-Poise) assuming an average van der Waals separation, 
h V j lv , of 3.4 A, the NT length, L m , of 53.65 A and the NT radius, R m , of 6.8 A. This result is close to the 
viscosity of pentane at 0.22 cP at room temperature [12], These values complete the parameterization of 
the interfacial model (4.7), which approximates the force-velocity dependence during the entire NT pull- 
out process (Figure 5). 



FIGURE 5. Comparison between the viscoplastic interface model and MD simulation data. 


6. Concluding Remarks. Carbon nanotube (NT) pull-out from a crystalline polymer matrix is 
simulated with molecular dynamics (MD) and an interfacial friction model. The resulting interfacial 
sliding between the adjacent surfaces of the NT and the surrounding polymer has been investigated with 
MD simulations and modeled by a nano-scale analog of Newton's friction law. Verification has been 
given that the length scales involved in the simulation ensure a meaningful representation of sliding 
properties. The friction model developed involves tire pull-out force, an effective viscosity, and the strain 
rate, which are parameterized directly from the MD data. 
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The entire pull-out process has been divided into four stages: 1) an initial equilibrium state, 2) the 
onset of NT slippage at the threshold pull-out force, 3) transition between the initial and steady sliding, 
and 4) steady sliding. When steady sliding is reached, a carbon-ring-based period in NT sliding velocity 
is observed, as the NT passes through the polymer matrix. On average, interfacial interactions during NT 
sliding can be approximated by a linear force-velocity relationship whose slope is the effective viscosity 
coefficient in the friction model. The magnitude of the viscosity is dependent on the chemical 
interactions at the NT/polymer interface, and is therefore, a potentially important structure/property 
relation for the design and optimization of novel NT-reinforced polymer composites. In this study, the 
weak interactions observed are inherent to the non-bonding van der Waals forces used to describe the 
NT/polymer interface. 

The interfacial friction model developed here is relevant to the evaluation of interfacial stress transfer 
in nanocomposites. 
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